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In this paper, we study the effects of attractive correlation on the topological Insulator {TI) with 
topological flat-bands using an extended attractive Kane-Mele-Hubbard model (KMHM). In the 
KMHM, we found a quantum phase transition from TI to the superconductor (SC) state upon the 
increasing of the attractive Hubbard interaction U at the mean field level. This type of SC phase 
transition is different from the traditional SC phase transition which develops from the gapless 
Fermi Liquid. Cooperon-type gapped excitations exist in the TI side near this type of SC phase 
transition. 


I. INTRODUCTION 

The integer quantum Hall (IQH) effect was first ob¬ 
served in a two dimensional (2D) electron gas subjected 
to strong perpendicular magnetic field [I|. This effect pro¬ 
vides the first example of topological states that beyond 
the Landau symmetry breaking paradigm. After this ob¬ 
servation, Haldane, in 1988, proposed a model (the Hal¬ 
dane model) and he found a state which also has the 
IQH effect in this model but the realization of this state 
doesn’t need the external magnetic fieldQ. The Haldane 
model describe a system spinless fermions and the time 
reversal symmetry is broken in this model by its complex 
NNN hoppings. The state that Haldane found is also a 
topological one. The two type topological states we men¬ 
tioned above that support IQH effect could be charac¬ 
terized by an topological invariant - TKNN number (the 
Chern number) [3j . Following the lesson from the Haldane 
model, people wonder naturally that whether the frac¬ 
tional quantum Hall (FQH) effect could also be realized 
in a model without external magnetic field. Recently, 
models with topologically nontrivial flat-bands (TFBs) 
were found to be an promising candidate to realize FQH 
effect without external magnetic field[^-Q. 

Along with the IQH and FQH states in which time- 
reversal symmetry breaking is required, time-reversal 
symmetry protected topological states of matters are also 
discovered in the quantized spin Hall effect (QSH)@, Q. 
People called them Z 2 topological insulator {TI) state. 
The typical model for the Z 2 topological insulator is the 
Kane-Mele (KM) modelQ. Recently, the correlated ef¬ 
fects uyTI states are studied by various groups with the 
Kane-Mele-Hubbard model as the starting point [ill [1^. 

In this paper, we investigate the effects of attractive 
interaction to a TFBs system using an attractive Kane- 
Mele-Hubbard model on the honeycomb lattice: FTkmh- 
There are two important parts in iJKMH : L^ekm and 
Hu, see Eq. ©■ In the original KM model, the authors 
generalizes Haldane’s model 0 to include spin with time 
reversal invariant spin-orbit interactions So the KM 
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model is a free model. FTekm is the free limit {U = 0) 
of iLxMH and is a general case of KM model i.e., the 
next-nearest-neighbor (NNN) hoppings for the spin-f and 
spin-j, electrons are complex valued and complex conju¬ 
gate to each other in iLEKM- So the Kane-Mele model is a 
special case of iLEKM i-e., when the complex valued NNN 
hoppings is purely imaginary or the Haldane phase (f>ij in 
Eq.(I2]) is 7r/2. There are particle-hole symmetry in the 
half-filing orignal Kane-Mele model i.e., when the com¬ 
plex valued NNN hoppings is purely imaginary it reduce 
to a spin-orbit interactions Q. So there are no particle- 
hole symmetry in iLEKM except 4>ij = 7r/2. If we varying 
the Haldane phase tpij there exist a so called TFBs limit 
in iJEKM [i-la, see FigHKb). 

Generally for the system with flat-bands, the kinetic 
energy will be quite suppressed and the interaction be¬ 
comes highly relevant. In this paper, Hu in Hkmh is 
treated by self consisted mean field method, at this level 
we find a, TI — SC quantum phase transition as the 
Hubbard interaction strength increases beyond a critical 
value Uc in the i/KMH model, see FigHl In this type SC 
phase transition, there is a kind of Cooperon-type exci¬ 
tations in the insulator side near the SC phase transition 
and the Cooperoii-tvpe excitations are gapped before its 
condensation HTHSI. 


II. TOPOLOGICAL FLAT-BAND MODEL 

The Hamiltonian that we study in this paper is iJKMH, 
the attractive Kane-Mele-Hubbard model which can be 
writed as: 

L^kmh = Hbkm + Hu — T ^ cl^Ca, (1) 

1,(7 

where a =t, 4- denotes the spin degree freedom, Hu = 
—U RitRil represent the on-site Hubbard type attrac¬ 
tive interaction, fj, is the chemical potential. We call 
Hekm the extended KM model whose Hamiltonian is 
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written as: 

HEKM = -t ^ ^ ^ (2) 

-1" Y. + ^■^■ 

{{{i,j))).o' 


In i^EKM, the first term describes the nearest-neighbor 
(NN) hopping on the honeycomb lattice and its hop¬ 
ping strength t is set as the unit of energy in the rest of 
this paper. The second term describes the next-nearest- 
neighbor (NNN) hopping with a Haldane type complex 
strength , the phase factor in this term is spin de- 

pendended, i.e., for spin f electron that hop¬ 

ping clockwise in the fundamental plaquette as shown by 
the blue arrow in Fig. 1(a) and for spin I elec¬ 

tron. This term recover the original spin-orbit coupling 
term of KM model when cj } = tt/2 (so here t' cannot be 
regarded as the the spin-orbit interaction strength). This 
term also reduces i?EKM’s full spin rotational SU(2) sym¬ 
metry to a U(l) symmetry. Thus there is time reversal 
symmetry (TRS) but no full spin rotation symmetry in 
.f^EKM- The last term represents the next-next-nearest- 
neighbor (NNNN) hopping with strength t", the energy 
bands of Hekm could achieve a flat-bands limit with this 
term in it (more about this limit see below). In this pa¬ 
per, we consider the half-filling case of Hkmh he., the 
chemical potential fi = 0. 

The Hekm model is the free limit {U = 0) 

of the Hkmh model. In terms of the basis vec¬ 


tor $'1 (fc) = 


?! /)f 


^fKMH could be ex¬ 


pressed in a block-diagonal matrix form as: Hekm = 
J2k (^) ^0 {k) $ (fc) where 


Hoik)=IoC + i\'^ (3) 

with h.)- = = (h“, —h®^), Iq is the 

4x4 unit matrice and cr is the Pauli matrice that 
act on the A, B sublattice space of the bipartite hon¬ 
eycomb lattice. The three components of vector 
= (Rey (fc), — Imy (fc), £>), where 7 (fc) = ^ (fc) -|- ^ {k) 
with C (fe) = -t Ei=i and C" {k) = -t" , 

withai = (l/2,-V3/2),a2 = (l/2, ^3/2), ag = (-1,0). 
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C = —2t cos (pij cos[/c • bi] 

i=l 

3 

D = —2t s'm (pij Y^ sin[fc • bi] (4) 

i=l 

with h = (3/2,-V3/2), 62 = (0,V3), 63 = 

(-3/2,-V3/2), see FigHKa). We set the lattice constant 
a = 1 in the rest of the paper. The energy dispersion of 
the Hekm can be founded by diagonalize Hq (fc): 


En{k) = C±^J (h-f + {hvf + {h^f. (5) 




FIG. 1: (color online) (a) The honeycomb lattice consists of 
two sublattices: A (the red sites) and B (the black sites). 
The NN (red), NNN (blue) and NNNN (green) displacement 
vectors are denoted by a, h and c respectively. The arrow in b 
also denote the positive phase direction of the complex factor 
that company the NNN hopping, (b) The TFBs in the 
energy spectrum of the Hekm with parameters: t jt = 0.6, 
t /t = —0.58, <\>ij = O.dTT. 


There are two energy bands in En{k), both of them are 
doubly degenerate due to the spin degree freedom and 
have a flat-bands limit, i.e., with parameters: (j) = 0.47r, 

t = 0.6t, t = —0.58t In this limit, the lower 

bands of En{k) will become flat, see FiglTJb). The flat¬ 
ness ratio of the flat-bands in this limit (the ratio of the 
band gap over bandwidth) can reach about 50. we know 
that the NN and NNNN hopping energy ^(fc) and ^"{k) 
are vanish at the two Dirac points: (27r/3, 27r/3-\/3) and 
(27r/3, —27r/3-\/3) in momentum space, but the NNN hop¬ 
ping energy (fc) dose not, just like the spin-orbit cou¬ 
pling term in KM model, it will opens an energy gap at 
those two Dirac points. So Hkmh is gapped in its free 
limit. We denote this energy gap as At', it’s the bulk 
gap of the TI state since before Hkmh model entering 
the SC phase it is in a T/ phase, the magnitude of Ap 
is related to t . Thus we called Af the bulk gap which 
playing the same role as the gap of semiconductor in la¬ 
in Hkmh’s Hatband limit, the Chern number for the 
spin up and down electron that filled the two lowest 
bands of i?„(fc) can be calculate in this way: 


Cf/t - 


47r 


n • 


/ dn dn \ 

\dkx ^ dky) 


cPk 


The Chern number of spin f components electron is = 
1 with n= h.^/|h.|.| and Cj, = —1 with n= h^/|h^|, 
reflecting the TRS in Hekm- The spin Chern number 
Cs = (C-f- — C\) /2 = 1 reflect the QSH in the free limit 
of Hkmh model [Q. 

This conclusion could be further verihed by the pres¬ 
ence of edge states in the energy spectrum of Hekm when 
it’s imposed a cylinder boundary condition with zigzag 
edges, i.e., we set periodic boundary condition along the 
system’s ^-direction and open boundary condition along 
the y-direction. The numerical results are depicted in 
FiglU from which we can see that there are topologi¬ 
cally protected edge states, which is one of the signa¬ 
tures of topological states of matter. So in half-filling 
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FIG. 2: (color online) The edge states of Hbkm in its flat- 
bands limit in a cylinder geometry with a zigzag edge. This 
is a typical edge states of TRS topological insulator which 
also reflecting the nontrivial topology of the flat-bands. 
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III. INTERACTING FLAT-BAND MODEL 


In this section, we study the effects of attractive Hub¬ 
bard type interaction of Hkmh rnodel in the flat-bands 
limit of iLfiKM with the mean field method. In this limit, 
the resulting mean field phase diagram of the ILkmh 
model could be see from FigUKa). From Fig0Ka) we could 
see that the TI phase in the Hkmh’s free limit is unstable 
against an SC phase transition as the interaction strength 
increase beyond a finite critical value: Uc — 3.03t. 

According to the self-consistent mean field method, 
firstly we introduce an SC order parameter fields to de¬ 
couple the Hubbard term Hu in Hkmh- The SC or¬ 
der parameter fields are defined as: As = (ci-|-Ci4,). By 
the mean field approximation, we substitute = 

Ci-i-Cij, = + S into the Hubbard 

term Hu in Eq.dl]), with ~ 

6 = {ci^Cii — {ci-fCii)) two small quantities. We discard 
the second order terms of and 6, then Hu could be 
decoupled into a bilinear form of and Cia as: 


Hu 


- H ^ — 

i 

+ Asc|'^cJj_ 


+ C/As|As|^ (6) 


The honeycomb lattice’s bipartite and translational in¬ 
variant, so we could introduce the usual Fourier trans¬ 
formation to the electron creation (destruction) operator 
cL (Go-) in F^kmh and we denote the fourier transform 
of 4GA,a and cleB,a as and bl^ respectively: 


FIG. 3: (color oline) The contour plots of the particle density 
distribution of the wave function of the zero-energy modes in 
coordinate space, this pair of zero-energy modes are bound to 
the TT-fluxes. The values of particle density are indicated by 
color in this figure,the parameters are in Hkmh’s flat-bands 
limit. 
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and flat-bands limit, with the two lowest flat-bands are 
filled, Hekm support the so called TFBs. Thus in the 
flat-bands limit of Hekm model, it is a Z 2 TI with TFBs 
as a consequence of the TRS in Hekm- 

On the other hand, the response of T/ to topolog¬ 
ical defects such as: 7r-fluxes has also been suggested 
as a probe of the nontrivial topology of the Z 2 TI 
M- Through numerical calculation, we really found that 
there are zero-energy modes in the energy spectrum if a 
pair of TT-fluxes are adiabatically inserted into the Hekm 
model in its flat-bands limit. From the particle density 
distribution in FiglSj we can see that those zero-energy 
modes are really located around the 7r-fluxes in coordi¬ 
nate space. Those results imply that the Hekm model’s 
topological non trivial in its flat-bands limit. 


where Ns is the number of unite cells and k belong to 
the first Brillouin Zone of the honeycomb lattice. Now 
we substitute Eq. dH]) and Eq. ([7]) into Eq. ([T]) we could got 
the momentum space form of Hkmh: Hkmh)^)- In the 
Nambu basis: ip^k) = {al^,a-ki,bl^,b-ki), Hkmh)^) 
could be casted into a matrix form in momentum space, 
that’s: 

HKMH{k)=Y,^Hk)h{k)i;{k) + UNs\As\^ ( 8 ) 

k 


with the 4x4 matrix h (fc): 


h{k) 


(C + D -UAs j{k) 0 \ 

-UAs -C-D 0 -7(fc) 

7t(/fc) 0 C-D -UAs 

\ 0 -"fHk) -UAs -C + D) 


(9) 
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FIG. 4: (color online) (a) The mean field results of As we solve 
from Eq. dlB at JTkmh’s flat-bands limit, (b) The excitation 
energy gap of the SC phase as a function of the interaction 
strength U, this excitation gap increases with U monotoni- 
cally so there is no further topological phase transition upon 
varying U. 



where C and D are the same as in Eq.Q. The quasi¬ 
particles spectrum can be found by diagonalize the 
Hamiltonian(|ni) in the momentum space: 

Er. (k) = ±Vei±2e2. (10) 

where 

e, = C^+D^ + \j{k)f + iUAsf, 
e2 = ^C2 (^D^ + \^{k)\^y 

We would investigate the instability of the TI phase of 
^^KMH in the presence SC fluctuation as the interaction 
strength U increases by minimizing the ground state’s 
energy Eq = J^nk^riik) against the SC order param¬ 
eters As, i.e., dEg/dAs = 0. Then the self-consistent 
mean field equation is given by: 

where Ei{k) = —\/ei -I- 2e2 and E 2 {k) = —^/ei — 2e2 are 
the two lowest filled bands in Ea. ifTOl) . 

The mean field solution of the self-consistent equation 
Ea. lfTlT) are plotted in Eigl4l(a) from which we can see 
that beyond a critical value of the attractive interaction 
strength: Uc — 3.03t the TI phase of the i/KMH model in 
its flat-bands limit become unstable against a SC phase 
transition. 

The quasi-particles excitation gap in the SC phase in¬ 
creases with U monotonically, see Eigl4l(b). This mean 
that there is no gap closing and no further topological 
phase transition upon varying U . 

FigE] and Figl6] is the mean filed results of the SC 
order parameter Ag we obtained from solving Eg. dill) and 
the excitation gap of SC phase over a range of t and 
U with the other parameters are fixed in the flat-bands 
limit. Note that in the presence of the bulk gap At' 
the order parameter and the excitation gap of SC 
phase are not the same like in the traditional SC phase 
transition[I3. From FigElwe could infer that the critical 


FIG. 5: (color online) The mean field results of the SC order 
parameter As we solve from Eq. GD) over a range of U and t , 
Other parameters are in the flat-bands limit. 

interaction strength Uc which separate the TI and SC 
phase increases nearly linearly with t . 

As a consequence of TRS in i/KMH, the Chern number 
will be zero and we verified this by a numerically calcu¬ 
lation. We calculate the Chern number of the two lowest 
bands oih(k) is /i(fc)[il[3 , the numerical results show 
that the total Chern number C = 0 in the SC phase 
atU = 3.03 with Ag = 0.0147 while the other parameters 
in the flat-bands limit. This mean that the SC phase 
transition destroy the TFBs. The topological property 
of this SC state can also be see from its edge excitations, 
so we also calculate the energy bands of T^kmh in a cylin¬ 
der geometry. The numerical results are depicted in Figl3 
from which we can see that the edge states are gapped 
when the SC order is developed in T/rmh- So beyond the 
critical interaction strength Uc, Hkmk entering a topo¬ 
logical trivial SC phase. On the other hand, we found 
a pair finite-energy bound state in the energy spectrum 
if a pair of 7r-fluxes are inserted into the SC phase with 
periodic boundary condition in both x and y direction, 
see FigEl comparing with the 7r-flux-induced zero-energy 
modes in Figl3]in the TI phase, show that the SC phase 
is topological trivial. This finite energy bound states and 
the gapped edge states in Fig [ 7 ] verified that the Hkmh 
model will loss its topology properties upon entering this 
s-wave SC phase as the attractive interaction strength 
increasing. 

So the SC phase in T/kmu develops from a fully gapped 
TI state. This kind of SC phase transition is different 
from the traditional SC phase transition which devel¬ 
ops from a Fermi liquid in which the low energy excita¬ 
tions are gapless, li a, SC order is develops from an fully 
gapped state, in our case an TI state, there will be an 
kind of gapped excitations, the so called Cooperon, in the 
insulator side of this SC phase transition[l3, [l^- This 
type of SC phase transition’s mechanism colud be under¬ 
stand in this way. There is a competition between the SC 
pairing gap Ag and the topological energy gap Af upon 
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FIG. 6: (color online) The mean field results of the excitation 
energy gap in the SC phase over a range of U and t while 
other parameters are in the fiat-bands limit. 


the varying of interaction strength U. If the topological 
energy gap A^/ is large enough, the Cooperon excitations 
will have an energy gap and will not condense, so the SC 
order parameter is vanishingly small. On the other hand, 
if Af is very small or the SC energy gap is large enough 
(at large t/), the Cooperon excitations will become gap¬ 
less and condensed, their condensation will leads to the 
SC phase transition at last. 




FIG. 7: (color online) Energy bands of Hkmh in its SC phase 
in a cylinder geometry with “zigzag” boundary. The SC 
phase is at [/ = 3.91 with As = O.Olfit, the other parame¬ 
ters are in the fiat-bands limit. 
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In this paper, we study the influence of attractive 
correlation on TI using the T/rmh model via a self- 
consistent mean field method. In the mean filed level, 
we found a TI — SC phase transition upon increasing 
the attractive interaction strength U in this model. The 
mechanism leading to this SC phase transition is differ¬ 
ent from the traditional ones. There is a competition be¬ 
tween the topological energy gap A^/ and the SC energy 
gap before the SC phase transition. This competition is 
reflect in the gapped Cooperon excitation, when the SC 
pairing gap is larger than A^/, the Cooperon excitation 
will condensed and the SC phase transition occurs in 
the system eventually. We think those results may help 
in better understanding the quantum exotic states in the 
correlated topological insulator and the pseudogap state 
in the high Tq cuprate suDerconductor(l^ . 


FIG. 8: (color online) The contour plots of the particle density 
distribution of the finite-energy bound states in the coordinate 
space. We can see that the finite-energy states are bound to 
the TT-fiuxes. The values of particle density are indicated by 
color. JTkmh’s in its SC phase and flat-bands limit. 
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